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Abstract 

The charmless bottom meson decays are systematically investigated based on an approximate six quark operator 
effective Hamiltonian from perturbative QCD. It is shown that within this framework the naive QCD factor- 
ization method provides a simple way to evaluate the hadronic matrix elements of two body mesonic decays. 
The singularities caused by on mass-shell quark propagator and gluon exchanging interaction are appropriately 
treated. Such a simple framework allows us to make theoretical predictions for the decay amplitudes with rea- 
sonable input parameters. The resulting theoretical predictions for all the branching ratios and CP asymmetries 
in the charmless B^, B^ , Bg tttt, ttK, KK decays are found to be consistent with the current experimental 
data except for a few decay modes. The observed large branching ratio in i? — > tt^tt" decay remains a puzzle 
though the predicted branching ratio may be significantly improved by considering the large vertex corrections 
in the effective Wilson coefficients. More precise measurements of charmless bottom meson decays, especially 
on CP-violations in B KK and Bg — > tttt, ttK, KK decay modes, will provide a useful test and guide us to a 
better understanding on perturbative and nonperturbative QCD. 



PACS numbers: 13.25.Hw,12.38.Bx,12.38.Lg,11.30.Er 



I. INTRODUCTION 



Hadronic B-meson decays play importance role not only for understanding the dynamical scheme of 
hadronic decays and testing the flavor structure of the Standard Model(SM), but also for probing the 
origin of CP violation and new physics signals beyond the SM. In particular, the precise measurement 
and systematic study for hadronic charmless B decays may provide a window for such purposes. The 
branching ratios of B ^ tttt and ttK modes have been measured with a good accuracy[l] and a 
large direct CP violation has been established in tt'^K' mode [1]. The most severe discrepancies 
between the experimental data and theoretical predictions come from the unexpected large branch 
ratio of B ^ tt^tt^ and some unclear CP violations in B ^ tt^K decays, which are called vrvr, ttK 
puzzles [2, 3]. Theoretically, to predict consistently those decays, it needs to deal with the short- 
distance contributions in a complete and systematic way from the high energy scale to a proper 
low energy scale at which the perturbative calculations remain reliable, and treat the long-distance 
contributions which contain the non-perturbative strong interactions involved in those decays. The 
main task is to reliably compute the hadronic matrix elements between the initial and final hadron 
states. Several novel methods based on the naive factorization approach (FA) and four quark operator 
effective Hamiltonian have been developed to evaluate the hadronic matrix elements, such as the QCD 
factorization approach (QCDF)[4], the perturbation QCD method (pQCD) [5], and the soft-collinear 
effective theory (SCET)[6]. These methods have been widely used in analyzing hadronic B-meson 
decays and made great progresses in understanding the hadronic structure and properties of strong 
interactions. To understand the puzzles whether they are due to the unknown new physics or it is 
because of the lack of our knowledge on the hadronic properties of strong interactions, it still needs 
to investigate further the various approaches within the framework of QCD and to check the validity 
of assumptions and approximations made in the practical calculations. 

The widely used theoretical framework of weak decays is based on the current-current four fermion 
operator effective Hamiltonian derived via operator product expansion and renormalization group 
evolution. In hadronic weak decays, the short-distance contributions of QCD are characterized by 
the Wilson coefficient functions of four quark operators and the long-distance contributions are in 
principle obtained by evaluating the hadronic matrix elements of four quark operators. The Wilson 
coefficient functions are in general calculated by perturbative QCD which is well developed, while 
the evaluation of hadronic matrix elements remains a hard task as it involves non-perturbative effects 
of QCD. To deepen our insights into the hadronic decays, we shall first reinvestigate the four quark 
operator effective Hamiltonian whether it is always suitable as a basic framework for all hadronic 
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weak decays. In fact, for the mesonic two body decays of B meson, it concerns three quark-antiquark 
pairs once each meson is regarded as the quark-antiquark bound state at the quark level structure. 
This fact then naturally motivates us to consider six-quark operator effective Hamiltonian instead of 
four-quark operator effective Hamiltonian. Namely, we shall begin with six quark diagrams of weak 
decays with both W-boson exchange and gluon exchange, and derive formally the six-quark operator 
effective Hamiltonian based on operator product expansion and renormalization group evolution when 
including loop corrections of six quark diagrams. We shall show how this approach allows us to figure 
out what are the assumptions and approximations made in effective four quark operator approach, 
and how the simple QCD factorization scheme can reliably be applied to evaluate the hadronic matrix 
elements with the six quark operator effective Hamiltonian. For the infrared singularity caused by 
the gluon exchanging interaction when evaluating the hadronic matrix elements of effective six quark 
operators, it is shown to be simply treated by the introduction of a mass scale motivated from the 
gauge invariant loop regularization method [7], where the energy scale fj,g is introduced to play the 
role of infrared cut-off energy scale without violating gauge invariance. 

The paper is organized as follows. In section II, after briefly reviewing the four quark operator 
effective Hamiltonian, we begin with the primary six quark diagrams with a single W-boson exchange 
and a single gluon exchange, and the corresponding initial six-quark operator. It is shown that a 
complete six quark operator effective Hamiltonian is in general necessary to include all contributions 
from both perturbative and non-perturbative QCD corrections, especially the non-pertubative QCD 
corrections at low energy scale /i < nic ~ 1.5 GeV could be sizable. To demonstrate how the six quark 
operator effective Hamiltonian provides a reliable framework for hadronic two body decays of B meson, 
we will focus, as a good approximation, on the dominant QCD loop diagrams of six quarks so as to 
avoid the tedious calculations. In section III, it is demonstrated how the QCD factorization approach 
becomes a simple and natural tool to evaluate the hadronic matrix elements of mesonic two body decays 
based on the six quark operator effective Hamiltonian. In particular, the so-called factorizable and non- 
factorizable, emission and annihilation diagram contributions are automatically the consequences of 
QCD factorization for the hadronic matrix elements of effective six quark operators. The treatment on 
the singularities caused by the gluon exchanging interactions and the on mass-shell fermion propagator 
is presented in Section IV. In Section V, all the amplitudes of charmless bottom meson decays are 
completely obtained by using the QCD factorization approach based on the approximate six quark 
operator effective Hamiltonian. Our numerical results with appropriate input parameters are presented 
in section VI, as a good approximation, the resulting predictions on branching ratios and CP violations 
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of charmless bottom meson decays are much improved and also more closed to the current experimental 
data. The conclusions and remarks are given in last section. The detailed calculations involved in the 
evaluation of various decay amplitudes are presented in the Appendix. 



II. EFFECTIVE HAMILTONIAN OF SIX QUARK OPERATORS 
A. Four Quark Operator Effective Hamiltonian 

Let us start from the four-quark effective operators in the effective weak Hamiltonian. The initial 
four quark operator due to weak interaction via W-boson exchange is given as follows for B decays 

Oi = {qibi)v-AiqjUj)v-A, = u, c, q'^ = d, s (1) 

The complete set of four quark operators are obtained from QCD and QED corrections which con- 
tain the gluon exchange diagrams, strong penguin diagrams and electroweak penguin diagrams. The 
resulting effective Hamiltonian (for b ^ s transition) with four quark operators is known to be as 
follows 

10 



^ q=u,c 



Cl{^l)o[''\^^) + C2{fi)oi''\f,) + c.(/x)o,(^^) 

i=3 



+ h.c. , (2) 



with Ag = VgbV*g and Vij the CKM matrix elements, Cj(/x) the Wilson coefficient functions[8] and 
Oi{^) the four quark operators 

O^l^ = {qibi)v~A{sjqj)v-A , O^^^ = {sibi)v--A{qjqj)v-A , 

O3 = {sibi)v-Aj2%^'j^v-^ ' ^4 = X(g-6i)v-A(sigj)y-A , 



(3) 



O5 = {sibi)v-A'Y{qjqj)v+A , Oe = -2'Y{qibi)s-p{sjqj)s+p , 

q' q' 

O7 = -^isibi)v-A'Yeg>{q'jqj)v+A , Og = -2>Yeq:{q[bi)s-p{sjq'j)s+p , 

3 ^' 3 ^' 

O9 = -j^isibi)v-AYeqi{qjqj)v~A , Oio = - ^ e<y'(g^6i)y_yi(sjgj)y-A , 

Here the Fermi constant Gp = 1.16639 x 10~^ GeV^^, and the color indices i, j, and the notations 
{q'q')v±A = q'7iii^^l5)q' ■ The index q' in the summation of the above operators runs through u, d, s, 
c, and b. The effective Hamiltonian for the b ^ d transition can be obtained by changing s into d in 
Eqs. (2)and (3). 
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B. Six Quark Diagrams and Effective Operators 

As mesons are regarded as quark and anti-quark bound states, the mesonic two body decays 
actually involve three quark-antiquark pairs. It is then natural to consider the six quark Feynman 
diagrams which lead to three effective currents of quark-antiquark. The initial six quark diagrams 
of weak decays contain one W-boson exchange and one gluon exchange, thus there are four different 
diagrams as the gluon exchange interaction can occur for each of four quarks in the W-boson exchange 
diagram, see Fig. 1. 




FIG. 1: Four different six quark diagrams with a single W-boson exchange and a single gluon exchange 

The resulting initial effective operators contain four terms corresponding to the four diagrams, 
respectively. In a good approximation, the four quarks via W-boson exchange can be regarded as 
a local four quark interaction at the energy scale much below the W-boson mass, while two QCD 
vertexes due to gluon exchange are at the independent space-time points, the resulting effective six 
quark operators are hence in general nonlocal. The six-quark operators corresponding to the four 
diagrams in Fig. 1 are found to be 

Og) = 4vra. //^^e-((--->+(---)'=)(,'(x3)7.TV(x3))- ' 



(2^)4 (2^)4- ^--^^'-^ ^ ^-^"k^ + ie 

(g-2(xi)ri 7"T"gi(x2)) * feCxOFagaC^i)), 



^2 ^ 7 7 (27r)4 (2^)4 - V. v-^.//.^ v^.;;^2+^g 

{q2{x2) / , . 7'-r"Figi(xi)) * {qA{xi)T2q3{xi)), 
p'^ — m^^ + le 



(27r)4 (2^)4- ^-^^^'^^ ^ ^^''fc2 + ie 

fe(xi)Figi(xi)) * (g-4(xi)F2 /^^'l- 7"T°g3(x2)), 



5^ ^ 7 7 (27r)4 (2^)4 ^ V. v-.//.^ v^.;y^2+^g 

'2 2 I ■ 

p^ — + ^£ 



fe(xi)Figi(xi)) * (g-4(x2) /^T"' .y T-T^q^ix^)), (4) 



where k and p correspond to the momenta of gluon and quark in their propagators, qi is usually 
set to be heavy quark like b quark, xi, X2 and X3 are space-time points corresponding to three 
vertexes. The color index is summed between qi,q2 and ^3,^4. Note that all the six quark operators 
are proportional to the QCD coupling constant due to gluon exchange. Thus the initial six quark 
operator is given by summing over the above four operators 

o(^)=x:o£). (5) 
3=1 

Actually, the initial six quark operators {j = 1,2,3,4) can be obtained from the following 
initial four quark operator via a single gluon exchange 

O = {q2^iqi) * iqiT2q3)- (6) 



C. Six Quark Operator Effective Hamiltonian via Perturbative QCD 

Based on the above considerations with the introduction of six quark operators, in this section 
we shall specify the initial six quark operator O^^^^^ {q = u, c) to the case of nonleptonic bottom 
hadron decays and show how to obtain six quark operator effective Hamiltonian. The initial six quark 
operator in b-decay with AS 7^ is as follows (for the b ^ d transition with AS* = 0, just replacing s 
hyd) 



1=1 

{ fe(xi)7^1 - 7^) ,^^'^']. YTtM^2m{x,h,{l - l')qj{x,)) 
— + le 



p'^ — + le 



+(g-.(xi)7^(l-7^)^.(xi))(.-,(x2) 7'^^fe.7M(l-7^)g.(^i))} 

p'^ — + le ■' 

k2 + ^^ (^'mixshuT^nlnixs)), (7) 
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which can be regarded as an effective operator resulting from the corresponding initial four-quark 
operator with a single gluon exchange 

O^i^ = {qibi)v-A{sjqj)v-A 

= ^7^(1 - l'')k){sj^,{l - j')q,) (8) 

with q = u, c. 

Similar to the procedure of obtaining the four quark operator effective Hamiltonian from the initial 
four quark operator O^^ of weak interaction, one should evaluate the six quark operator effective 
Hamiltonian from the initial six quark operator O^^^^'' when running the energy scale from m\Y to the 
low energy scale fi ~ mb- As the first step for finding out the complete set of independent effective six 
quark operators, one needs to evaluate all possible one loop diagrams based on the initial six quark 
diagrams (Fig. 1). The possible six quark diagrams at one loop level are plotted in Fig. 2. 

It is useful to classify those diagrams into three types: type I is the loop diagrams in which only 
the effective four quark vertex of weak interaction receives loop corrections including the penguin 
type loops, the single gluon exchanging interaction for six quark operators remains mediating between 
one of four external quark lines of loops and a spectator quark line (see Fig. 2a); type II is the 
loop diagrams where only the single gluon exchanging vertexes receive loop corrections (see Fig. 2b); 
the remaining loop diagrams are regarded as type III in which one of the gluon exchanging vertexes 
touches to the internal quark/gluon line of loops (see Fig. 2c). Note that in Fig. 2a and Fig. 2b we 
only plot, for an illustration, the six quark diagrams with a gluon exchanging between one of the four 
external quark lines of effective weak vertex and a spectator quark line, while for each of them, there 
are actually three additional different diagrams corresponding to other three choices of external quark 
lines, they are omitted just for simplicity. 

To evaluate all the diagrams is a hard task, as a good approximation, we shall pay attention to 
the type I and type II diagrams. The type III diagrams are in general suppressed at the perturbative 
region with energy scale around mf, as they involve more internal quark lines and contain no large 
logarithmic enhancements. From the evaluation of four quark operator effective Hamiltonian, it is 
known that when the energy scale runs via the renormalization group evolution from the high energy 
scale at /i ~ mw to the low energy scale around /i ~ mi,, the loop corrections of type I diagrams 
should result in the six quark operators with all effective four quark operators and the corresponding 
Wilson coefficient functions, meanwhile the loop corrections of type II diagrams will lead the strong 
coupling constant of the gluon exchanging interaction to run from high energy scale at mw to the 
low energy scale at /u. Thus, when ignoring the type III diagrams, we arrive at an approximate six 
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FIG. 2: The diagrams in (a) arc loop contributions only to the effective weak vertex (type I), and diagrams in 
(b) are loop contributions only to the giuon vertexes (type II). The diagrams in (c) are loop contributions for 
both weak and strong vertexes (type III). 

quark operator effective Hamiltonian as follows 

V ^ j=l q=u,c 
10 

+ Y,>^f^CdM%{^^)} + h.c. + ..., (9) 

with the CKM factor Xg^'^^ = ^b^*^^^)- The dots represent other possible terms that have been 
neglected in our present considerations. O- (j = 1,2,3,4) are six quark operators which may 
effectively be obtained from the corresponding four quark operators Oj(/u) (in Eq. (3)) at the scale 
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fi via the effective gluon exchanging interactions between one of the external quark hnes of four 
quark operators and a spectator quark hne at the same scale ^. The general forms and definitions 
of ^i^q (a*) U — 1)2,3,4) for the corresponding four quark operators Oi{fi) are similar to the ones of 
Oq^^ {j = 1,2,3,4) given in Eq. (4) but with replacing as(miy) by Usin) due to QCD corrections of 
type II diagrams. 

Before proceeding, we would like to point out that a complete six quark operator effective Hamil- 
tonian may involve more effective operators from the type III diagrams and lead to a non-negligible 
contribution to hadronic B meson decays when evaluating the hadronic matrix elements of six quark 
operator effective Hamiltonian around the energy scale fi ~ y^2AQCD^b ~ "Ttc ~ 1-5 GeV where the 
nonperturbative effects may play the role. We shall keep this in mind and regard the above six quark 
operator effective Hamiltonian as an approximate one. 



III. QCD FACTORIZATION BASED ON EFFECTIVE SIX QUARK OPERATORS 

We shall apply the above effective Hamiltonian with six quark operators to the nonleptonic two 
body decays of bottom mesons. The evaluation of hadronic matrix elements is the most hard task 
in the calculations of the decay amplitudes. In this section, we are going to demonstrate how the 
factorization approach naturally works for evaluating the hadronic matrix elements of nonleptonic 
two body decays of B meson with six quark operators. 

To be explicit, we here examine the hadronic matrix element of -B ^ vr'^vr" decay for a typical six 
quark operator oj^^^ 

(6) ^ f f (fk (fp .((^^„^^)p+(^3_^3)fc) 1 1 

J J (27r)4 (27r)4 - mj 

which is actually a part of the six quark operator O^^^ in the effective Hamiltonian. Its hadronic 
matrix element for B tt^tt^ decay leads to the following most general terms in the QCD factorization 
approach 



--< vrVO 



oS I Bo > 



(^k (^p ^-i((xi-X2)p+(x2-X3)fc) 1 



(27r)4 (27r)4 



p2 



0^0 



< TT vr 



[dk{x2){i> + md)YT^il^{l - 7')&i(^l)][di(xi)7^(l - -f'')d,{xi)][dra{x3huT^„dn{x3)] \ Bq > 



.0(1) 



.0(2) 



rO(3) 



rO(4) 



(11) 



with 



/1^C)(2) ^ f f " ^ »P -,-^((xl-X2)p+(x2-X3)fc) 1 ^g^g 



f (fk 




(27r)4 


(27r)4 


7^(1 - 




f d^k 




(27r)4 


(2^)4 


Y(i - 




f (fk 




(2vr)4 


(2vr)4 


'7^'(1 - 




( d^k 






(2^)4 


V(l - 


7')]p. 



^jO(3) ^ / / "'^ "'^ g-»((^i-X2)p+(x2-X3)fc) 1 ^gyg 



dl 



7i^O(4) ^ /" f ^-»((xi-X2)p+(x2-a:3)fc) 1 ^g^g 

LL II /9^H ('9^^4 k'^{p'^-ml) " 



3^,^(xi,X2)M^f:^(xi,X3)Mi;j( 



where M-^^^{xi,Xj) = [Mx{xi,Xj)]nm {X = B,7r) with n, m the color indices and a, (3 the spinor 
indices, is the hadronic matrix element of two quark operators for a single meson X. In light-cone 
QCD approach, it is found to be [9] 



M/„Z^{x„x,) = <n<^[P)\d'^^{x,)di{x,)\0>='-^^-^ l\xe-^(^^^^+^^^^^^ (13) 

with Fm {M = B, tt) the decay constants. Here M^°'{u, Pb) and M^"(x, P) are the spin structures for 
the bottom meson and light meson vr and characterized by the corresponding distribution amplitudes 

M^"(n,PB) = -[mB + rBf(I^B{u)]pa, 
M^-{x,P) = [rr'<t>Ax) - f^nlH€i^) - io^,rrf<g{x) + ia^.P^ '^''^^^ QkJ '^^''' ^^^^ 

with V = ^-j*^ 1 n = n'^ + n~ — V and (fP" = (j)"' /Q. The light-cone distribution amplitudes (pfjiu) 
{M = B,7r, X = —,p,T) are given in [9] up to twist-3. The definition of momentum for quarks and 
mesons is explicitly shown in Fig. 5. As a good approximation, both the light quarks and light mesons 
are taken to be massless, i.e., P^ = 0. 

It is interesting to note that the four amplitudes M^j^^ {i = 1,2,3,4) are corresponding to four 
diagrams (l)-(4) in Fig. 3. The first diagram is known as the factorizable one, the second is the non- 
factorizable one and color suppressed. The third is the factorizable annihilation diagram and color 
suppressed, and the fourth is an annihilation diagram and its matrix element vanishes. 
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FIG. 3: Different ways of reducing hadronic matrix element of effective six quark operator by QCD factorization 
approach. 



After performing the integration over space-time and momentum, the above amphtude is simphfied 
to be 

1 .1 .1 1 



with 



M^LiBTTTT) =< 7r°7r° \ofl\Bo>= [ [ [ du dx dy 

Jo Jo Jo 



(uP*-(l-x)Pi) 



2 



m 



d 



(Pi -uP+Y - ml ((1 - x)Pi + yP2-uP^) 
{xPi+P2y-my {xPi + {l-y)P2-uP+y-my' ^ ^ 



Tr[M^(y,P2)7'^(l-7')] =^^i^B F^Mn)mlf,My)€{^), 



((1 - x)f>i+ y f2- u md)7^(l - 7^)] 
^'"^Fb Fl(l)B{u)m%{mB{u + x + y- 2)(j)^{x) + ^^(1 - x){4P^{x) - 4>l{x)))(j)^{y), 



^lI = F^2T¥[M,(x,Pi)7.M^(y,P2)7"(a;n+f2+m,)7;.(l-7')]Tr[MB(t/,PB)7^(l-7')] 

= '^Fb F^MuWB{xmlMy)M^) + 2^' ((1 + ^)€{^) - (1 - x)4>l{x)W^{y), ) 
M^^l = 0*Fb F^Tr[M^{x, Pih,M^{y, P2h^{l - 7') 

MB{u,PB){xfi+ (1 - y)r2-ulfpr mrf)7V(l - 7')] = 0, (16) 

where M^j^ {i = 1,2,3,4) are obtained by performing the trace of matrices and determined by the 
distribution amplitudes. Cp = resulted from summing over the color indices. It can be seen 

that M^j^ corresponding to Fig. 3.(1) is color allowed, M^j^ and M^j^ corresponding to Fig. 3. (2) and 
Fig. 3. (3) are color suppressed, while M^j^ corresponding to Fig. 3. (4) vanishes as it is not allowed for 
colorless mesons. 
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From the above explicit demonstration, it can be seen that the simple QCD factorization approach 
becomes a natural tool to evaluate the hadronic matrix element of effective six quark operators in the 
mesonic two body decays. For a given effective six quark operator, its hadronic matrix element for 
mesonic two body decays gets four different combinations in the QCD factorization approach, namely 
it consists of four different amplitudes corresponding to four topologically different diagrams. From 
the above example, it is noticed that the amplitude M^^^ is a color-allowed factorizable one in an 
emission diagram, M^j^^ is a color-suppressed non-factorizable one in an emission diagram, M^^^ 
is a color-suppressed factorizable one in an annihilation diagram, while m'^^^ vanishes as it cannot 
match to a colorless meson. 

When generalizing the above analysis to the present framework based on the approximate six quark 
operator effective Hamiltonian, there are in general four types of six quark diagrams corresponding 
to four types of effective six quark operators, their hadronic matrix elements for two body mesonic 
decays lead to sixteen kinds of diagrams (see Y\gA.{ai)-[di), i = 1,2,3,4) as each of the effective six 
quark operators leads to four kinds of amplitudes in the QCD factorization approach. 




FIG. 4: Four types of effective six quark diagrams lead to sixteen diagrams for hadronic two body decays of 
heavy meson via QCD factorization. 

It is known that the effective four quark vertexes concern three types of current-current interactions: 
{y -A)x{V - A) or (LL), {V - A) x {V + A) or {LR), (5 - P) x (5 + P) or thus each of the 
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diagrams in Fig. 4 actually contains three kinds of diagrams corresponding to three types of current- 
current interactions. Therefore, there are totally 48 kinds of hadronic matrix elements involved in the 
QCD factorization approach, while it is easy to check that only half of them are independent with the 
following relations: 









— tFt in ■ 




- A^r IN ■ 




= 0- 




^LRa^ 


^^^LR 




^"LR 


= A^r, IN^- 


^"LR 


= 0; 


Mgp 


rpF . 
— ^SPa^ 


Mgp 
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]\/ra3 

Mgp 
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Mgp 


n. 
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_ rpF . 
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Mil 


= nLb/Nc, 


Mil 


= ^iLb/Nc, 


Mil 


= 0; 




_ rpF . 


Mil 


= TfpJN,; 


Mil 


= ^SPJ^C, 


Mil 


= 0; 


Mfp 




Mfp 


= Tf!nb/Nc; 


Mfp 


= ^iRb/Nc, 


Mfp 


= 0; 


Mil 


= 0; 


Mil 


= T^UNc, 


Mil 


= ^IlJNc; 


Ml\ 




Min 


= 0; 


Mfn 


= TfpjN,; 


Mil 


= ^SPa/Nc, 


Mil 


— ^LRa-- 


Mfp 


= 0; 


Mfp 


= npjN,; 


Mfp 


= ^IrJNc: 


Mfp 


- A^ ■ 


Mil 


= 0; 


Mil 


= TlJN,; 


Mil 


= A'lLb/Nc; 


Mil 


- A^ ■ 


Mfn 


= 0; 


Mil 


= TfpjN,- 


Mil 


= AsPb/Nc, 


Mil 


- A^ ■ 




= 0; 


Mfp 


= T[nb/Nc; 


Mil 


= ^^Rb/Nc; 


Mil 


- A^ 



(17) 



where T^^ and T^j^ {X = LL, LR, SP) represent the factorizable emission diagram contributions, 
and T^j^ {X = LL, LR, SP) are the non-factorizable emission diagram contributions. A^^, A^^ 
and A^^, A^^ {X = LL, LR, SP) denote the so-called factorizable and non-factorizable annihilation 
diagram contributions respectively. Their detailed definitions and general formalisms are presented in 
the Appendix. 



IV. TREATMENT OF SINGULARITIES 



In the evaluation of hadronic matrix elements, there are two kinds of singularities, one is caused 
by the infrared divergence of gluon exchanging interaction, and the other arises from the on-mass 
shell divergence of internal quark propagator. As the quark propagator singularity is a physical-region 
singularity, one can simply add ie to the denominator of quark propagator and apply the Cutkosky 
rule [10] to avoid such a singularity. It then allows us to obtain the virtual part of amplitudes as 
the Cutkosky rule gives a compact expression for the discontinuity across the cut arising from a 
physical-region singularity. In general, a Feynman diagram will yield an imaginary part for the decay 
amplitudes when the virtual particles in the diagram become on mass-shell, and the resulting diagram 
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can be considered as a genuine physical process. It is well-known that when applying the Cutkosky 
rule to deal with a physical-region singularity of all propagators, the following formula holds 

1 



— ml + ie 



— m'^ + ie 



P 
P 



p2 — rn^ 



-iTrS[p'^ - ml], 



p"^ — Trig 



-m<5[p2-m2], (18) 



which is known as the principal integration method. Where the first integration with the notation of 
capital letter P is the so-called principal integration. 

For the infrared divergence of gluon exchanging interactions, only adding ie to the gluon propagator 
is not enough as such an infrared divergence is not a physical-region singularity, one cannot simply 
apply the Cutkosky rule. To regulate such an infrared divergence, we may apply the prescription used 
in the symmetry-preserving loop regularization[7] which allows us to introduce an intrinsic energy scale 
without destroying the non-abelian gauge invariance and translational invariance. The description of 
the loop regularization is simple: evaluating the Feynman integrals to an irreducible integrals, replacing 
the integration variable and integration measure J by the regularized ones via [7] 



with conditions 

N 



lim^^cf'(M;2)" = 0, c^ = l (i = 0,l,--- ,iV and n = 0,l,---), (20) 

where cj^ are the coefficients determined by the above conditions. With a simple form for the regulator 
masses Mi = fig + I Mr (/ = 0, 1, • • • ), the coefficients is found to be = (— 1)^ {n~i)\ i\ ' ^^sA, 

which leads the regularized integrals to be independent of the regulators. Here the energy scale 
Mq = fig plays the role of infrared cut-off but preserving gauge symmetry and translational symmetry 
of original theory. 

In the present case, there is no ultraviolet divergence for the integral over k as it is constrained by 
the finite momentum of hadrons, so all the terms with / 7^ in the summation over I vanish in the 
limit M/j — > 00. As a consequence, it is equivalent to add an intrinsic regulator energy scale iig in the 
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denominator k'^ in Eq. (11), thus one can use the usual principal integration method to avoid such a 
singularity, i.e., 

1 1 „r 1 1 . ...o 2 



k"^ k"^ — + ie 



i7r6[k^ - fi^]. (22) 



With the above considerations, the singularities appearing in the integrations over k and p can simply 
be avoided by the following prescription 

11 1 1 



(23) 



^2 ^p2 _ j^2~) ^^2 _ ^2 _j_ ^p2 _ j^2 _|_ ' 

Note that as the gauge depending team k^k^ can be transformed to the momentum |^ on exterior 
line of spectator quark, they are all on mass shell in our present consideration (as defined in Fig. 5 in 
the appendix), their contributions equal to zero, thus our results are gauge independent. 
Applying this prescription to the amplitude illustrated in previous section, we have 

.1 .1 rl 1 



m2l{B^t^) =< ^°7r° I O'fl \Bo>= [ [ [ dudxdy 

Jo Jo Jo {U^B 



{uP+ -{l-x)Pi)^-fil + ie 



(Pi -uP^Y -ml^ie {[l-x)Pi+yP2-uP^Y -ml + ie 

_l hk I LL f24^ 

{XP1 + P2Y -ml + ie {xPi + {l-y)P2-uP^Y -ml + ie^' ^ ' 

V. AMPLITUDES OF CHARMLESS BOTTOM MESON DECAYS 

With the above considerations and analyses, the QCD factorization approach enables us to evaluate 
all the hadronic matrix elements of nonleptonic two body decays of B meson based on the approxi- 
mate six quark operator effective Hamiltonian. The amplitudes of charmless B meson decays can be 
expressed as follows: 

A{B^ ^ ^+^-) = V,,,V:,\Pr{B) + \p'i^{B) + PF(i?) + 2PY{B) + \p^^ {B) - ^A^-^{B)] 

-VudV:,[T^^B) + E^^B)], 
A{B+ ^ vr+vrO) = -L{y^,y;^[p|^(s) + pf-(i?)] - V^dV:,[T-^ [B) + C^-{B)]}, 

A{B^^7r\^) = -L{-y,,y;jp-(i?)-P-^(i?)-ipf-(i?) + P-(S) + 2Pr(5) 

+ip#^-(P) - ^pE^^B)] + VudV:,[-C^'^iB) + E-^B)]}, (25) 
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for B — > vrvr decay amplitudes, and 

^(iJ+^vr+^O) = -VtsV,l[P¥^{B)-^PE{^^{B) + P^^{B) + ^PE{f{B)] + VMA-^\B), 

A{B+^7r'K+) = l={VM[P¥''{B) + P^^{B) + ^PE^^{B)+P^^<{B) + ^Pi^^{B)] 

-VusV:,[T^''{B) + C'^'^iB) + A-^(S)]}, 
A{B'^7r-K+) = VM[P¥''{B) + ^PE^''{B)+PE^iB)-^P§^^{B)]-VusV:,T-^{B), 
A{B'^^'K') = -l^{VuVl[P^^{B)-P^^{B)-\p%^{B) + Pl^{B) 

-\pE^^{B)] + VM^'^'iB)}, (26) 
for B — > -kK decay amplitudes, and 

A{B' ^ K+K-) = -VM * [P^iB) + Pf ^(S) + \P^§^{B) - + VudV:,E^^{B), 

A{B+^K+K') = -VM[P^^iB)-^PE{^^B)+P^^{B) + ^PE^'<m+VudV:,A^^{B), 
AiB^^K^K'') = -VuV:,[P^^{B)-^P%^{B)+pf''{B) + P^^{B)+P^''{B) 

-IpZ^'iB) - \p^^Hb) - \pEVm, (27) 
for B KK decay amplitudes, and 

^(S°^7r+7r-) = -VME^^B,) + VM[2PriBs) + lpi^^B,)], 
^(^O-vrV) = i=A(i?o^vr+vr"), 

^(5° - n+K-) = VM[P§-{B,) + ^Pf|-(5.) + Pf -(5.) - ^^11^(5.)] - KsK*63^^"(^s), 
^ 7r°K°) = --L{y,,y;,[pf -(i?,) - Pf^(i?,) - ^Pf^-(i?.) + P|-(i?.) - IpE^'^iBs)] 
+VMC^^{Bs)], 

A{B^,^K+K-) = -Vt,V:,[P§^{Bs) + \p^^^{B,) + P§'<{B,)+P^^{B,) + P^^{B,) 

+lpi§''{B,) - Ipi^^'iB,) - \pE^'HBs)] + VusV:,[t'^'^{B,) + E^^iB,)], 
A{B',^K'K') = -VM[P§''{B,)-\pE^^{Bs) + P§^{B,)+P^^{B,) + P^^{B,) 

-IP^^^'iBs) - \pi§^{Bs) - \pE^^{B,)1 (28) 



for Bs vrvr, vrK, KK decay amplitudes. The eleven types of amplitudes T^'^'^^{M), 
C^hM2^M), P^^'^^^M), P^^^^M), A^-^i^'^2(M), ^A'^iA^2(Af), P^^^^'^^M), P^^'^^^M), P^^'^^^M), 
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p|^i*^2(Af), P^^^^^iM), with M1M2 = TnT,TTK,KTT,KK,KK are defined as follows 



for the so-called emission diag 
A^^^^{M) = 47ra,(/i 



'NM1M2 

LL 



(M) 



+ 

ams, and 
Gp 



(29) 



NM1M2 

SP 



pEmM2^M^ = 47ra,(/. 

+ [C8(/^ 

for the so-called annihilation diagrams. Where T^j^, '^xAi ^XA^ ^XA = LL,LR,SP, A = a,b) 
arise from the hadronic matrix elements and their detailed expressions are given in Appendix. Note 
that T^^^ and T^^^ are slightly different as the wave functions of K meson and K meson are not 
equal at high order in the twist expansion. 



NM1M2 

LL 



(M) 



^-^..,'^LR W}> (30) 
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When redefining the above amplitudes to the widely used diagrammatic amplitudes in the phe- 
nomenological analysis, 

T = T''''{B), C = C'^'iB), E = E''''{B), P = + PW{B), Pa = ^P^ , 

Pew = PEwiB)^ Pew — Pew (-^)' Pew — Pew iP)^ Pew = Pew{P)^ 

T' = T^^{B), C = C^^{B), A' = A^'^iB), P' = P:}^^{B) + P^^{B), P'^ = 2Pf ^, 

P'ew — PEwiP)i P'ew — Pew^ {P)i P'ew — Pew^{P)^ P'ew — Pew^{P)^ (^1) 
P" = P^^iB) + P^^iB), P'i = P^^iB) + P^iB), P'^ = PE§^iB), A" = A^^iB), 
Pe^t = [Pi^^iB) + Pi^^{B)]/2, P^i, = [Pi^^iB) - P#|^(B)]/2, P^^ = pEV{B). 

the decay amplitudes can be reexpressed in terms of the familiar forms in the diagrammatic decom- 
position approach 

VuVliP + Pa + \p^w + \pew - \pEw) - VudV:,{T + E), 

^[Vuv:,{Pew + pEw) - yndv:,{T + c)], 

-^[VuV,U-P -Pa + Pew + \pEw + \pew " \pEw) " K^KT^CC - E)], 
VtsVliP' - \p'Ew + \p'iw) + ^usK^A!. 

^ytsViip' + \p'Ew + p'ew + \pew) + yusv:,ij' + c + ^o], 

1 / 1 / I ^ IE I 

-^j^tsVt\kP - -^Pew ~ Pew ~ -^Pew) + ^usV*hC ], 
VMiP'i + ^P^i, + P^i.) + VudV:,E", 

vMiP" - \p'k^ + Ip'^^) + v^,v:,^', 

VMiP" - ^P'jfy, +P'i- ^P^i. - ^P^U (32) 

It is noticed that there is a slight difference to the usual diagrammatic decomposition approach with 
the extra contributions from the annihilation electro- weak diagrammatic amplitudes P^^^^^^ and 
Pew^^^^ ^ which are actually small and neglected in the usual diagrammatic decomposition approach. 



A{B^ - 


-> vr^vr ) 


A{B+ ■ 


vr+vrO) 


A{B^ 




A{B+ - 


vr+^O) 


A{B+ - 




A{B° - 




A{B° - 


ttO^O) 


A{B^^ 


K^K~) 


A{B+ - 




A{B° - 
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A similar redefinition can be made for Bg decays, 



Eg — E'"'^{Bs),PsA — '^PA'iBs),PfEW — Piwi^s), 

n = T^^B,), C = C^^B,),P', = P§^{B) + P§^Bs),P'sEw = P§wiBs), 

^sEW — ^EW \^s),rg^^r — f^iY [-Ds), 

T': = T^^iB,), P'J = Pi^iB,) + P§^{B,), P'J^ = P^^iB,) + P^^(Bs), P'JEw = PE^^'iB), 
E': = E^^iBs),P'Ji^ = [Pi^^iBs) + Pi|^(i?.)]/2, P'Ji^ = [P^^^iB,) - Pi^^{Bs)]/2, 
P'sEW = Pew^{Bs). (33) 

Then the decay amphtudes can be reexpressed as follows 

^(i?,^7r+^-) = vuv:,{PsA + \p^Ew)-yudv:,E,, 

A{Bs ^ ttV) = -^AiBs ^ vr+TT-; 
v2 

A{Bs 



v2 

a{Bs^t:-k+) = vuv:,{p', + '^pi%w-\p'J'Ew)+yusv:x. 

A{B, ^ vrOi^O) = -ljy,,Vl{Pl - \PI%^ - P^EW - IPs^i^ 
A{B,^K+K~) = VuVUP'jA + ^P'Jiw + P'Jiw) + VudV:,El 

A{B, K'K') = vuv:,{p': - \p':Ew + p'Ia - \p':iw - \p's 



-^P'sEW ~ -^P'sEw)- {'^'^ 



VI. NUMERICAL CALCULATIONS 



We are now in the position to make numerical calculations. 



A. Theoretical Input Parameters 

The short distance contributions characterized by the Wilson coefficient functions for the effective 
four quark operators were calculated by several groups at the leading order(LO) and next-to-leading 
order(NLO) [11], their values mainly depend on the choice for the running scale ^. In our numerical 
calculations, it is taken to be 

PL = ^2kQCDm^ ~ (1.5 ± 0.1)GeV. (35) 
The Qs value in the six quark operator effective Hamiltonian is also taken at /i = (1.5 it 0.1) GeV. 
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When considering the NLO Wilson coefficient functions and as, one needs to include the magnetic 
penguin-like operator Osg which is defined as [8] 

Osg = ^mb%a^,(l + 75)I^^G'^^''6, , (36) 

where i, j are the color indices. The magnetic-penguin contribution to the B ttK, tttt decays leads 
to the modification for the Wilson coefficients corresponding to the penguin operators, 



a^M '^ifiW--^ yr^Csgin) (37) 



with Cg^ = Cgg + C5 and |/^| = Where 04^6 are known to be defined as 04^6 = C^^q + which 

appear in the factorizable diagrams. 

For other parameters, we take the following typical values 

ruB = 5.28GeV, = 139.6MeV, m^o = 135MeV, rrib = 4.4GeV, = 1.5GeV,m, = O.lGeV, 

mu = md = 5MeV,/ij = 216.19MeV, U = ISO.lMeV, Fk = 159.8MeV,;U^ ~ 1.7GeV, 

UK - 1.8GeV,TB0 = 1.536ps, tb+ = lM8ps, A = 0.2272, A = 0.806, p = 0.195, f] = 0.326. (38) 

Especially, for the infrared energy scale fig introduced in this paper to regulate the infrared divergence 
of gluon exchanging interactions, we take the typical value of fig to be a universal one around the 
hadronic bounding energy scale of non-perturbative QCD 

fig = (400 ± 50)MeV. (39) 

To evaluate numerically the hadronic matrix elements of effective six quark operators based on the 
QCD factorization, it needs to know the twist wave functions of mesons. For the wave function of B 
meson, we take the following form [12] in our numerical calculations: 

2" 



4)b{x) = Nbx (1 — x) exp 



1 / xme Y 

'2\ujB J 



(40) 



For the light meson wave functions, it needs to know the twist distribution amplitudes which 
contains the twist-2 pion (kaon) distribution amplitude (l)Tt{K)^ aiid the twist-3 ones 4'^(^k-j and (fj^^^K)^ 
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they are parameterized as [13]: 




4>7t{K){x) 





1)^ -30(2j; - 1)2+3) 



(41) 



(42) 



(1 - lOx + lOx^) 



(43) 



In our numerical calculations, the shape parameters in the distribution amplitudes are taken the 
following typical values: 



iUB = 0.25 GeV ,u;b, = 0.33 GeV , rjs = 0.015 , = -3 , 

al = , af = 0.06 , af = 0.10 ± 0.10 , = 0.15 ± 0.15 , of =< = ± 0.10 (44) 



where the shape parameters for the bottom mesons are taken from [12], and other shape parameters 
are taken to fit the data. Since those shape parameters can vary by 100%, they agree with the ones 
in Refs.[13, 14], All parameters for the light mesons are taken at the energy scale 1 GeV [15], run to 
our calculation scale. Note that they may vary significantly when the scale runs to different values. 

B. Numerical Results and Discussions 

The numerical results for the CP averaged branching ratios and CP violations of charmless B meson 
decays are presented in Table I for B vrvr, ttK decay channels and in Table II for B KK channels. 
In Table III, we give the results for the branching ratios and CP violations for Bg vrvr, ttK, KK decay 
channels. The LO and NLO are corresponding to the leading order hadronic matrix elements with 
the leading order and next-to- leading order Wilson coefficients (which include the magnetic penguin- 
like operator Osg)- For comparison, the numerical predictions from the QCDF approach and pQCD 
approach are also listed in the Tables. It is seen that most resulting predictions in our present 
calculations are in good agreement with experimental data within the possible uncertainties from 
both experiments and theories, while it remains unclear how to understand the puzzles in the decay 
channel B^ — > vr'^vr'^ for its large branching ratio and possible positive CP violation, and in the decay 
channel B 7r^K~^ for the unexpected large positive CP violation. 
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As shown in[13], adding vertex corrections may improve the CP violation. The vertex corrections 
[16] were proposed to improve the scale dependence of Wilson coefficient functions. The ref.[16] 
has considered vertex corrections that only influence the Wilson coefficients of factorizable emission 
amplitudes. Those coefficients are always combined as C2n~i + %f and C2n + and modified by 



4tt 



Nr 



Nr 



Nr 



n 



1 



(45) 



with M2 being the meson emitted from the weak vertex. In the NDR scheme, Vi{M) are given by [16] 

121n(^^) - 18 + /q dx(/>A/(x) g{x) , for i = 1 - 4, 9, 10 , 
Vi{M) = I -i2ln{^) + 6- f^dx4>M{x)g{l-x) , fori = 5, 7, (46) 
—6 + dx 4'\j{x) h{x) , for i = 6, 8 , 

'Pm{x) / (f^j^iix) is the twist-2/twist-3 meson distribution amplitudes defined in Eq. 14. The functions 
g{x)/h{x) used in the integration are: 



9ix) 



6 I — mx — ITT 



+ 



1 — X 



2Li2(x) - In^x + 



21nx 
1 - X 



(3 + 2i vr) In X — (x 1 — x) 



h(x) = 2 Li2(x) — In^ x — (1 + 2i vr) In x — (x <-> 1 — x) 



(47) 
(48) 



Such a correction does not include the contributions of the ffist two diagrams in Fig. 2a which are 
considered as a part of form factor or meson amplitude. It is interesting to notice that the vertex 
corrections do improve the predictions for CP violations and bring CP violations in the decay channels 
7r°7r° and B+ tt^K^ to be more close to the experimental data. 

To enlarge the branching ratio of i? ^ tt'^tt'^, we shall examine an interesting case that only two 
vertexes concerning the operators Oi and O2 receive additional large nonperturbative contributions, 
namely the Wilson coefficients ai = Ci + ^ and a2 = C2 + jf are modified to be the effective ones: 



a2 - = C2(/i) + ^(/^) + ^C'f^(^2(M2) + Vo) , 



(49) 



Taking the value Vq = 25, the resulting branching ratio for B tt^tt^ becomes consistent with the 
experimental data. 
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It is more interesting to consider the possible nonperturbative effects by taking tlie effective color 
number Nc^'^ in color-suppressed diagrams. The numerical results with Nc'^'^ = 1.7 are presented in 
Table IV- VI, which provides an alternative explanation to the puzzle of observed large branching ratio 
B TT^TT^. For comparison, we also list in Table IV- VI the predicted results via the S4 scenario in 
QCDF [17]. 

The method allows us to calculate the relevant transition form factors at maximal recoil (with NLO 
Wilson coefficients including magnetic penguin contribution), 

TpB^TT _ n ncn+0. 029+0.10 rpB-^K _ n 009+0.034+0.013 rpBs^K _ r, 07/1 +0.023+0.013 /c^ 
— U.^O^_o. 024-0.010' -^0 — "-'•"J^^-0.029-0.011 ' -^0 ~ ' ^-0.013+0.0005 ' W"-"^ 

with input parameters /ig=400MeV, /i=1.5GeV, ;U7r=1.7GeV, ;UA'=1.8GeV.The first error arises from 
the range for fig = 350 ~ 450 MeV, the second error is caused by the running scale /i = 1.4 ~ 1.6 
GeV. 

The resulting form factors agree with the ones obtained from the light-cone QCD sum rule of heavy 
quark effective field theory [18] 

F(f ^'^ = 0.285j:J];[!J^, F^~*^ = 0.345 ± 0.021, F^''^^ = 0.296 ± 0.018, 

and from the full QCD sum rule[19] 

F(f ^'^ = 0.258 ± 0.031, F^^^ = 0.331 ± 0.041. 

To know the relative contributions from various diagrams and hadronic matrix elements of effective 
six quark operators, we present in the Table VII and Table VIII the numerical results for different 
kinds of topology amplitudes, the predictions for the strong phases are all relative to the leading order 
tree amplitude phase 6t — 1-93 in B ^ vrvr decay. It is interesting to see that the amplitudes from 
the annihilation diagrams are significant in comparison with the color suppressed emission diagrams. 

The predictions of 5*^0 in each method are almost the same and obviously larger than averaged 
data in PDG. But some new data in [20, 21] prefer a larger prediction. 

VII. CONCLUSIONS 

Based on the approximate six quark operator effective Hamiltonian derived from perturbative QCD, 
the QCD factorization approach has naturally been applied to evaluate the hadronic matrix elements 
for charmless two body decays of bottom mesons. The resulting predictions for the decay amplitudes, 
branching ratios, and CP asymmetries in B^, B^, Bg — > vrvr, ttK, KK decay channels have been 
found to be consistent with the current experimental measurements except for a few decay modes. 



23 



TABLE I: CP averaged branching ratios and CP violations for B tttt, ttK decay channels. The central values 
are obtained with parameters: /ig=400MeV, /i=1.5GeV, fi^^lJGeV, pij^^l.SGcV. The first error arises from 
the range for fig — 350 ~ 450 MeV, the second error stems from the running scale fi ~ 1.4 ~ 1.6 GeV. 



Mode 


Data [1] 


QCDF[17] 


pQCD[13] 


This work 


LO NLO(+MP) 


LO NLO(+MP) 


B+ ^ TT+A-o 
B+ ^ n"K+ 
B" -> n-K+ 
B" ^ 7r"A'" 


23.1 ±1.0 
12.8 ±0.6 
19.4 ±0.6 
10.0 ±0.6 


10 Q+l. 9+11. 3+1. 9+13. 2 
^y-«J-1.9-7.8-2.1-5.6 

1-1 -1+1.8+5.8+0.9+6.9 
^-'^•-'^-1.7-4.0-1.0-3.0 

-ifj q+2. 6+9. 6+1. 4+11. 4 
-'^"■•->-2. 3-6. 5- 1.4-4.8 

7+0.7+4.7+0.7+5.4 
' -0.7-3.2-0.7-2.3 


17.0 24.1 
10.2 14.0 
14.2 20.5 
5.7 8.7 


16.46 21.60tl:itt36 

9.12 11.78t^;81+2;22 

14.42 19.03t^j°t|;^^ 
6.61 ^Mtlf,t\Z 


B° 7r-n+ 
B+ ^ tt+tt" 
B" n^n" 


5.16 ±0.22 
5.7 ±0.4 
1.31 ±0.21 


o Q+4. 0+3. 6+0. 6+1. 2 
°-^-3.4-3.0-1.0-0.8 

r. r,+3.0 + 2.1+1.0 + 0.4 

"•"-2.4-1.8-0.5-0.4 

n o+0. 2+0. 2+0. 3+0. 2 
'-'•'J-0. 2-0. 1-0. 1-0.1 


7.0 6.7 
3.5 4.1 
0.12 0.29 


6 63 c 71+1.69+0.70 

O.DO 0. ( i_j 24-0.57 

4.43 4.69lJ:??l"o:^2"6 

0.11 o.i6l°:rrol 


Acp{7T+K°) 

Acp{n-K+) 


0.009 ±0.025 
0.047 ±0.026 
-0.095 ±0.013 
-0.12±0.11 
0.58±0.17[20, 21] 


n nnQ+"-002+o. 003+0. 001+0. 006 

U.UUy_Q QQ3_Q QQ3_Q QQ]^_Q 

n 071 +"-01'i'+0-'320+0. 008+0. 090 
"•'-" ^-0.018-0.020-0.006-0.097 

n n/l'i+0-011+0.022+0. 005+0. 087 
U.U+0_o 011-0.025-0.006-0.095 

n n'?'?+0-010+0-013+0. 005+0. 034 
U.UOO_o 008-0.016-0.010-0.033 


-0.01 -0.01 
-0.08 -0.08 
-0.12 -0.10 
-0.02 0.00 
0.70 0.73 


0.016 +o.oi6;S:gg^;o".oT 
-0.093 -o.o8ot°:°°«t°:r3 

150 -0 194+0-014+0.008 
U.iOU U.lZ+_o. 014-0. 007 

-0 006 -0 noi-"-OOi-0-003 
u.uuo u.uui_|_o. 000+0. 002 

711 n 71 5-° °i2-o.oo3 

u./ii rd_|_o 002-1-0.003 


Acp{n^7:^) 
Acp{n+n") 
Acpin^n^) 


0.38 ±0.07 
0.04 ±0.05 
0.36 ±0.32 
-0.61 ±0.08 


n nfir^+0-021+0. 030+0. 001+0. 132 
U.UOd_Q 021-0.028-0.003-0.128 

n nnn+0-000+0. 000+0. 000+0. 00 
u.uuu_o. 000-0. 000-0. 000-0. 000 

n /I CI +0.184+0.151+0.043+0.465 
U.+Oi_o. 128-0. 138-0. 141-0. 616 


0.14 0.20 
0.00 0.00 
-0.04 -0.43 
-0.34 -0.41 


0.178 0.187t°:gS?t°o:°olt 

0.000 o.oool°:™n:ro 

571 -0 r;47+0.018+0.046 
U.O/i U.0+(_o 025-1-0.033 

528 -0 cifii -0.011-0.010 
u.ozo u-ooi^o. 011+0. 009 



TABLE IL B A' AT modes with the same input parameters as Table L 



Mode 


Data [1] 


QCDF[17] 


pQCD[22] 


This work 
LO NLO(+MP) 


B+ ^ if+i^o 
B° K+R- 


1.36 ±0.28 
0.96 ±0.20 
0.15 ±0.10 


1 oc+0. 45+0. 72+0. 14+0. 91 
^•'J'J-0. 39-0. 49-0. 15-0.40 

1 qr+0. 41+0. 71+0. 13+1. 09 
^•«J"J-0. 36-0. 48-0. 15-0.45 

n ni Q+0. 005+0. 008+0. 000+0. 087 
U.UiO_o. 005-0. 005-0. 000-0. Oil 


1.65 
1.75 


0.85 imtlf+_ViA 

0.65 oMtr^tr^ 

0.07 o.o7iro^lE5:"ol 


Acp{K+RO) 
Acp{K^R°) 
Acp{K+R-) 


0.12 ±0.17 
-0.58 ±0.7 


n -iGo+0. 047+0. 050+0. 016+0. 113 
U. 100^0.037-0. 057-0. 017-0. 133 

n i«7+0. 047+0. 045+0. 015+0. 046 
U.IU/ -0.037-0. 051-0. 017-0. 036 




0.096 o.o78tg:°}n:^°ol 

000 noo+"-""o+0-000 
u.uuu u.uuu_o. 000-0. 000 

n sn7 n o^o-o.ooo-o.ooo 

tJ^^^^^o. 042+0.000 
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TABLE III: Bg mr, ttK, KK modes with the same input parameters as Table I. 



Mode 



Data [1] 



QCDF[17] 



pQCD[14] 



This 



LO NLO(+MP) 



Bs 
Bs 
Bs 
Bs 
Bs 



t:+K- 

K+R- 



0.5 ±0.5 



5.0 ± 1.25 



24.4 ±4. 



0.024 
0.012 



-0.003+0.025+0.000+0.163 
-0.003-0.012-0.000-0.021 

+0.001+0.013+0.000+0.082 
-0.001-0.006-0.000-0.011 

ir, 9+4.5+3.8+0.7+0.8 
^'-'•^-3.9-3.2-1.2-0.7 

n /iQ+0. 28+0. 22+0. 40+0. 33 
'-'■^^-0. 24-0. 14-0. 14-0. 17 

99 7+3.5+12.7+2.0+24.1 
' -3.2-8.4-2.0-9.1 

r,A 7+2.5+13.7+2.6+25.6 
' -2.4-9.2-2.9-9.8 



0.57 
0.28 



+0.16+0.09+0.01 
0.13-0.10-0.00 

+0.08+0.04+0.01 



7.6 



-0.07-0.05-0.00 

+3.2+0.7+0.5 
-2.3-0.7-0.5 



n 1 c+0. 05+0. 10+0.02 
"•^"-0.04-0.05-0.01 



13.6 



+4.2+7.5+0.7 
-3.2-4.1-0.2 



1 ^ G+5. 0+8. 3+0.0 
-'^'^•"-3. 8-4. 7-0.0 



0.19 
0.10 
6.96 
0.07 
13.26 
15.25 



23+''-'^-'^+°-"''' 
0.11 



-0.01-0.05 

+0.01+0.03 
"^-0.01-0.02 



7 n9-l-l.ll-f0.63 
' •'-'^-0.91-0.51 



0.09 
16.68 



+0.04+0.03 



0.03-0.02 

+5.37+4.32 
-3.71-3.24 



IS? Q/t+5. 80+4.56 



Acp{tt+tt ) 

Acp{7:+K-) 

Acp{K+K-) 
Acp{K"K^) 



-0 oi2+'^-'^"-'-+°-'^-'-^+"''^"-^ 
-0.012 



-0.004-0.012-0.001 

+0.001+0.012+0.001 
-0.004-0.012-0.001 



0.39 ±0.17 



_n nfi7+0-021+0. 031+0. 002+0. 155 
U.UO I -0.022-0.029-0.004-0. 152 

n /I ifi+0. 166+0. 143+0.078+0.409 
"•*l"-0. 120-0. 133-0. 145-0.510 



r, 9^1+0.039+0.033+0.023 
^•^^1-0.036-0.030-0.012 



0.594 



-0.018+0.074+0.022 
-0.040-0.113-0.035 



0.040 
0.009: 



-0.010+0.020+0.005+0.104 
-0.010-0.023-0.005-0.113 



-0.002+0.002+0.001+0.002 
-0.002-0.002-0.001-0.003 



n 9"? Q+0. 009+0. 049+0. 008 

-U.ZO.O_g QQ2_g Q44_o oil 



0.018 

0.018 
0.182 
0.128 
-0.218 
0.000 



n m c+o. 028-0.003 

u.uitJ_Q Q20+0.002 



0.015 



+0.028-0.003 
-0.020+0.002 



n 1 oq+0. 012+0. 018 

U.100_Q ggg_o 015 

n nc;/i-t-" oi4-i-o.089 

^•"3^-0. 014-0. 081 

_n 1 or;+0.014+0.007 
U.iOO_g 010-0.009 

n nnn+0 000+0.000 

U.UUU_Q ggg_Q ggg 



The puzzles for the observed large branching ratio in — > decay and possible large positive CP 
violations in B ^ irK^ decay need to be further investigated. As we have emphasized at the beginning 
that the six quark operator effective Hamiltonian considered in this paper is an approximate one, and 
a large number of six quark diagrams which suppressed at high energy scales have been ingored, but 
they may become sizable at low energy scales. Furthermore, when given the predictions, we have 
only considered the uncertainties caused by the choices of running scale /i and infrared energy scale 
fig as their effects are more significant than others. In general, the theoretical uncertainties could be 
much larger when the possible uncertainties for all the input parameters are included. The masses of 
light mesons are also neglected in comparison with the bottom meson masses, i.e, m^/m^ and 
m\/m?^ 0. 

Nevertheless, it is remarkable that such a simple theoretical framework based on the approximate six 
quark operator effective Hamiltonian from the perturbative QCD and the naive QCD factorization for 
the nonperturbative QCD effects can result in a satisfactory theoretical prediction for the charmless 
B meson decays Bg — > tttTjTtK, KK . It also shows that the singularity due to the on mass- 
shell fermion propagator can simply be treated with the principal integration method by apply the 
Cutkosky rules, and the one caused by the gluon exchanging interactions can well be regulated by the 
description used in the loop regularization method with the introduction of an intrinsic energy scale 
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TABLE IV: The same as Tabic I but including the vertex contributions and compare with QCDF S4. 



Mode 


Data [1] 


QCDF S4[17] 


pQCD[13] 


This work 








LO 


NLO±Vertex 


LO 


NLO±Vertex 


eff 


c 


B+ ^ 


23.1 ± 1.0 


20.3 


17.0 


OA c + 13.6(+12.9) 
8.1 (- 7.8) 


16.45 


99 nR + 7.39+4.25 
ZZ.UU_4 gg_3 21 


22.06 


19.50 


B+ ^ tt"K+ 


12.8 ±0.6 


11.7 


10.2 


-,o Q+10-0(+ 7.0) 
5.6 (- 4.2) 


9.12 


-,9 r)ri+3-84+2.19 
lZ.UU_2.54_i.65 


11.66 


10.95 


B" ^ 7T-K+ 


19.4 ±0.6 


18.4 


14.2 


90 Q + 15.6(+11.0) 
8.3 (- 6.5) 


14.41 


1Q 09+6.67+3.84 
iy.JZ_4 4i_2.gi 


19.62 


18.68 


B" ^ 7r°K" 


10.0 ±0.6 


8.0 


5.7 


q ,+ 5.6 (+ 5.1) 
3.3(- 2.9) 


6.61 


a QC+3.25+1.88 
0.»0_2.14-1.42 


9.70 


8.71 


B° 7r-7r+ 


5.16 ±0.22 


5.2 


7.0 


a c+ 6.7 (+ 2.7) 
"•"J- 3.8(- 1.8) 


6.62 


7 n7+l-67+0.71 
' -1.29-0.58 


5.38 


4.89 


B+ ^ 7r+^" 


5.7 ±0.4 


5.1 


3.5 


4 0+ (+ 1-7) 
1.9(- 1.2) 


4.43 


A 97+0.96+0.33 
^■^ ' -0.73-0.29 


6.98 


6.43 




1.31 ± 0.21 


0.7 


0.12 


r, 9q+0.50(+0.13) 
"•^"^-0.20 (-0.08) 


0.11 


n 1 O+0.07+0.05 
"•^°-0.04-0.03 


1.03 


0.98 




0.009 ±0.025 


0.003 


-0.01 


-0.01 ±0.00 (±0.00) 


+0.016 


U.UZU^O. 003+0. 001 


0.020 


0.018 




0.047 ±0.026 


-0.036 


-0.08 


„ p., +0.03 (+0.03) 
"•"^-0.05 (-0.05) 


-0.093 


n n'^t^+o-006+0.004 

U.UOO_o. 002-0. 002 


-0.068 


-0.0529 




-0.095 ±0.013 


-0.041 


-0.12 


„ p,„+0.06(+0.04) 
"•""-0.08 (-0.06) 


-0.150 


n 1 oo+0. 015+0. 008 

U.100_Q gii_Q gQ7 


-0.117 


-0.131 




-0.12 ±0.11 


0.008 


-0.02 


„ p,_+0.03(+0.01) 
"•^ -0.03 (-0.01) 


-0.006 


n nc^i +0.003+0.000 

U.U01_o. 002-0. 000 


0.002 


-0.029 




0.58±0.17[20, 21] 




0.70 


PI -o+0.03(+0.01) 
' "^-0.02 (-0.01) 


0.710 


n 71 n-0. 002-0. 001 
U. ( 1U_Q QQ2^o.002 


0.789 


0.781 




0.38 ±0.07 


0.103 


0.14 


q ,„+0.20(+0.07) 
'^•^°-0.12 (-0.06) 


0.178 


n 1 oo+0. 002+0. 015 
U.100_o. 002-0. 014 


0.214 


0.223 


AcpilT+TT^) 


0.04 ±0.05 


-0.0002 


0.00 


0.00 ±0.00 (±0.00) 


0.000 


n nnn+0 000+0.000 

U.UUU_Q QQQ_Q QQQQ 


0.000 


0.000 




0.36 ±0.32 


-0.19 


-0.04 


n fiq+0.35(+0.09) 
"■""^-0.34 (-0.15) 


-0.571 


n /i7n+o oio+o o32 

"■'^"-'-0.011-0.018 


-0.174 


-0.208 




-0.61 ±0.08 




-0.34 


n .0+1 00 (+0.05) 

"•^"^-0.56 (-0.05) 


-0.528 


n cc,f;-0. 010-0.009 

U.OOU^Q QQ4^Q QQg 


-0.586 


-0.479 



TABLE V: The same as Table II but including the vertex contributions and compare with QCDF S4. 



Mode 


Data [1] 


QCDF S4[17] 


pQCD[22] 


This work 










LO 


NLO±Vertex 


eff j^eff 
"1,2 ^^c 


B+ ^ K+K° 


1.36 ±0.28 


1.46 


1.65 


0.85 


-1 1 O+0.26+0.18 
-L- J^'J-0.17-0.14 


1.13 0.85 




0.96 ±0.20 


1.58 


1.75 


0.65 


n Q7+0.22+0.16 
"■°'-0.14-0.11 


0.87 0.608 


B" K+K~ 


0.15 ±0.10 


0.070 




0.07 


n 07+0.03+0.01 
-0.10-0.01 


0.07 0.29 


Acp{K+K^) 


0.12 ±0.17 


-0.043 




0.096 o.o8ol°:"o;,^t"o:™? 


0.080 0.207 




-0.58 ±0.7 


-0.115 




0.000 o.oooiror^o-.r 


0.000 0.000 


Acp{K+K-) 








0.807 0.842;° °°^;°:^ 


0.84 0.78 



Hg. In particular, it is found that such a scale takes a typical value fig = (400 ± 50) MeV which is 
around the binding energy of hadron due to non-perturbative QCD effects. 

We would like to point out that although the theoretical framework discussed above is a much 
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TABLE VI: The same as Table III but including the vertex contributions and compare with QCDF S4. 



Mode 



Data [1] 



QCDF S4[17] 



pQCD[14] 



This work 



LO 



NLO+Vertex 



a 



eff 
1,2 



Bs - 


TT+TT 


Bs - 




Bs - 


-> 11+ K- 


Bs - 




Bs - 


-> K+R- 


Bs - 





0.5 ±0.5 



5.0 ± 1.25 



24.4 ±4. 



0.155 
0.078 
8.3 
0.61 
36.1 
38.3 



0.57 



+0.16+0.09+0.01 
-0.13-0.10-0.00 



n 9O+0-08+0.04+0.01 
"■^"-0.07-0.05-0.00 



7.6 



+3.2+0.7+0.5 
-2.3-0.7-0.5 



n 1 05+0. 10+0.02 
"■^"-0.04-0.05-0.01 

-, o f.+4.2+7.5+0.7 
-'^'-'■"-3.2-4.1-0.2 

1 «+5. 0+8.3+0.0 
-'"^■0-3. 8-4. 7-0.0 



0.19 

0.10 
6.96 
0.07 
13.26 
15.25 



0.23 



+0.01+0.07 
0.01-0.05 



n -I -r +0.01+0.03 
"■^^-0.01-0.02 

7 qr + 1. 15+0.63 
' ■•-'^-0.94-0.51 

r, -1 7+0.04+0.04 
"■^ ' -0.03-0.03 

1 R 77+5.36+4.23 
' ' -3.69-3.17 

18 04+5-72+4.34 
i0.a'±_3 gg_3 26 



0.23 

0.11 

5.73 

0.69 

16.97 

18.94 



0.69 
0.34 
6.58 
0.60 
15.76 
16.63 



AcpiTT+n ) 

Acp{tt+K-) 

AcpiK+K-) 
Acp{K"K") 



0.39 ±0.17 



0.109 
0.046 
-0.047 
0.006 



— n ni 9+0-001+0.012+0.001 

'-'■"^^-0.004-0.012-0.001 

n r,1 O+O-OOl+0. 012+0. 001 
U.Ui^_Q QQ4_Q gj^2-0.001 

241+*'-°^^+°-°'^^+*' °^'^ 



0.594 



-0.018+0.074+0.022 
-0.040-0.113-0.035 



n 9"^ Q+0. 009+0. 049+0. 008 

-U.ZO.O_g QQ2_g Q44_g oil 



0.018 

0.018 
0.182 
0.128 
-0.218 
0.000 



n m c+o. 028-0. 003 

U.UitJ_Q oj^g^o go2 



+0.028-0.003 
-0.019+0.002 

-0.009+0.015 
-0.002-0.016 

+0.017+0.017 
-0.011-0.006 



0.015 
0.182: 
0.831 



_n 1Q/1+0-014+0.010 

"•-'^^-0. 011-0. 010 

n nnn+o-ooo+o.ooo 

U.UUU_o ogo-O.OOO 



0.015 
0.015 
0.207 
-0.135 
-0.168 
0.000 



0.016 
0.016 
0.171 
0.057 
-0.191 
0.000 



TABLE VII: Diagrammatic amplitudes relating to the CKM matrix element A„ with 10 ^ GeV. 



topology 


T 


C 




A 


E 






LO 


81.48 


4.931e- 


-0.75i 




7.329e- 


3.07i 


TTTT 


NLO(+MP) 


81.46 


6.952e- 


-0.50i 




7.321e' 


3.07i 




NLO+Vertex 


83.51eO-04' 


13.88e- 


-1.60i 




7.321e- 


3.07i 




LO 


100.0 


6.020e^ 


-0.75i 


39.86e-"-50» 






ttK 


NLO(+MP) 


100.0 


8.558e- 


0.488i 


39.58e-''-5i* 








NLO+Vcrtcx 


102.4e"-04' 


17.52e- 


-1.59j 


39.58e-°-5°* 








LO 








2.267e2-3"* 


7.169e- 


0.85i 


KK 


NLO(+MP) 








2.045e2-3i* 


7.088e- 


0.85i 




NLO+Vcrtcx 








2.045e2-3i* 


7.088e- 


0.85i 



simplified one, it turns out that as the first order approximation the six quark operator effective 
Hamiltonian considered in this paper can be taken as a good starting point. We have actually examined 
two interesting cases by considering teh effective Wilson coefficient functions the effective color 

number Nc^^ in the color suppressed diagrams to bring the prediction for the branching ratio B —> 
be consistent with the experimental data. It is of interest to calculate high order contributions 
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TABLE VIII: Diagrammatic amplitudes relating to the CKM matrix element with 10 ^ GeV. 



topology 


Pt 




Pew 


pC 
^EW 


Pa 




pA 


pE 
^EW 


P = Pt + Pe 




LO 


6.555e- 


3.10i 


1.176e- 


-3.13i 


0.101e°-5°* 


2.210e°-°^* 


3.137ei-5i* 


0.021e- 


2.65i 


O.llle' 


-0.62i 


6.987e2-^2j 


TTTT 


NLO(+MP) 


7.478e- 


3. Hi 


1.175e- 


3.13i 


0.076eO-69' 


2 442e0 06» 


3.339ei-5i» 


0.019e- 


2.59i 


0.112e- 


0.63j: 


7.876e2-74i 




NLO+ Vertex 


7.684e- 


3.07i 


1.201e- 


3.09i 


0.232eO-93' 


2 442e0 06» 


3.339ei-5i» 


0.019e- 


2.59i 


0.112e- 


0.63i 


7.971e2'^8j 




LO 


8.332e- 


3.10i 


L490e- 


-3.12i 


0.123eO''9^ 




4.604ei-^** 






0.165e- 


-0.58i 




TlK 


NLO(+MP) 


9.483e" 


3. Hi 


1.489e- 


3.12i 


0.092e°-6S' 




5.069ei-9i' 






O.lSOe- 


o.eoi 


12.00e2'^6^ 




NLO+ Vertex 


9.731e- 


3.07i 


1.521e- 


3.07i 


0.287e°-96j 




5.069ei-9i* 






O.lSOe- 


0.60i 


12.11e2-79' 




LO 


10.57e- 


3.08i 






0.161e°"4* 


1.572e°-5°* 


3.433ei-69* 


0.069e- 


1.26i 


0.053e- 


^3.00i 


11.30e2-90' 


KK NLO(+MP) 


12.03e- 


3.08i 






0.120e°-59' 


1.667e0-54j 


3.663ei-^0j 


0.063e- 


1.31i 


0.052e- 


2.94i 


12.79e2-92, 




NLO+ Vertex 


12.35e- 


3.05i 






0.362eO-95' 


1.667e0-54' 


3.663ei-7''' 


0.063e- 


1.31i 


0.052e- 


2.941 


22.9ge2-95i 



TABLE IX: Diagrammatic amplitudes relating to the CKM matrix element A„ in Bs decays with 10 ^ GeV. 



topology 


Ts 


Cs 






Es 






LO 












2.28i 


TTTT 


NLO(+MP) 










3.980e- 


2.28i 




NLO+Vcrtex 










3.980e- 


3.28i 




LO 


77.75e"-0S' 


6.221e- 


0.84i 


53.96e2-58* 








NLO(+MP) 


77 74g0.08» 


7.876e-^ 


0.60i 


53.76e2-59* 








NLO+Vertex 


79.63e0"8'' 


14.96e- 


1.54i 










LO 


95.41e""**' 


8.085e- 


0.84i 


2.705ei-^" 


7.313e- 


1.44i 


KK 


NLO(+MP) 


95.41eO"8* 


lO.lOe- 


0.61i 


2.478ei-^^* 


7.219e- 


1.43i 




NLO+Vertex 


97.73eO"**'' 


19.36e- 


l.SOi 


2.478ei-^8» 


7.219e- 


1.43i 



though it is a challenging task. On the other hand, the precise measurements of charmless bottom 
meson decays, especially the measurements on CP-violations in B KK and Bg vrvr, ttK, KK 
decays, will provide a useful test for various theoretical frameworks. It is expected that more and 
more precise experimental data in the future super B-factory and LHCb will guide us to arrive at a 
better understanding on perturbative and nonperturbative QCD. 
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TABLE X: Diagrammatic amplitudes relating to the CKM matrix element At in Bg decays with 10 ^ GeV. 



topology 






PsEW 


pC 
^sEW 


PsA 








Ps = PsT + PsE 


LO 


- 




- 


- 


2.365ei™* 


- 


0.026e- 


2.24i 






2.265ei-"" 


TTTT NLO(+MP) 










2.597e°-99» 




0.025e- 


2.23i 






2.597eO-99^ 


NLO+ Vertex 










2.597e0-99» 




0.025e- 


2.23i 






2.597eO-99^ 


LO 


6.329e- 


3. oil 


1.122e-306i 


0.108e°-^9* 




4.105eii2» 






0.044e- 


2.78i 


5.315e2-5S' 


■kK NL0(+MP) 


7.2156" 


3. Oil 


1.121e-306j 


0.088ei"0' 




4.362ei05' 






0.037e- 


2.25i 


5.745e2-62i 


NLO+ Vertex 


lAlOe- 


3.02i 


1.145e~3.05j 


0.243ei-''^' 




6.275ei-22» 






0.093e- 


1.45i 


7.180e2-35' 


LO 


8.231e- 


3.00i 


0.918e+305i 


0.143e°-^^* 




5.852ei-62^ 


0.083e- 


1.47i 


0.043e^ 


1.63i 


9.671e2-63i 


KK NLO(+MP) 


9.348e- 


3.01i 


0.917e+3°5» 


0_ii8eO.95j 


1.320ei-5''j 


6.238ei-62j 


0.078e- 


1.48i 


0.047e- 


1.61i 




NLO+ Vertex 


9.597e- 


3.02i 


0.937e+304» 


0.307ei"4» 


1.320ei-5''' 


6.238ei-62' 


0.078e- 


1.48i 


0.047e- 


1.61i 


10.80e2-66i 
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APPENDIX: CALCULATIONS OF HADRONIC MATRIX ELEMENTS 

In this appendix, we are going to present the explicit expressions for all the hadronic matrix 
elements evaluated from the naive QCD factorization method based on effective six quark operators. 
To be specific, we shall first make the following convention for the momentums of quarks and mesons, 
which is explicitly shown in Fig. 5 

Where we have ignored the light quark mass in external lines and light meson mass to simplify 
calculation. 

Let us first give the factorizable emission contributions for the {V—A) x {V—A) and iV—A) x (y+A) 
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FIG. 5: Definition of momentum in B ^ MiM2. The ligfit-cone coordinate is adopted with (n"*", n ,k±) 



effective four quark vertexes, they are simply denoted by LL and LR 

-1 rl rl 



TlLa^'{M) = i-^^ Fm Fm.FmJ I I dudxdymUMiu) 

c i/ 

{mB{2mb - 7nBx)<pMAx) + ^lmA'^^bx - jn^) [(/>^^^ (x) - 0fj^(x)]}(/)M2(2/)^ra(^'2;), 
1 Cp 

'^LU^'^i^) = ^2-^ Fm^ j j j dudx dym%fiMi<pMiu)(j)M2iy)(pM^ix)hTbiu,x), 

r™i^^(M) = -r™^^^^(M). (A.i) 

The factorizable emission contributions for the {S — P) x (S + P) effective four quark vertex are found 
to be 

Fspa^'^i^) = '^t:^ Fm Fmi Fm2 [ [ [ dudxdy ms fJ-Ahff'Miu) 
^ ^^c Jo JO JO 

{mB{2mB - mb)(l)Miix) + fi]\i^[Anib - {x + (j;) + HM^msil - x)(j)lj^{x)} 

T^Pb'^'iM) = 4^ Fm Fm, Fm2 C C C dudxdy m| m^Mu) 

^ '^c Jo Jo Jo 

[niBUcpMj (x) + 2(1 - u)i^LM^4>M^ (^)]'/'a/2 iy)^Tb{u, x). (A.2) 
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Similarly, we obtain 



rpNMiM2 
rpNMlM2 

''^LLa 



'''LLb 



(M) 
(M) 

(M) 



"1 rl rl 

L ^ p ^ ^ / 

-I 



F f f f 

4 Jo Jo Jo 

{(n - y)mB(l)Md^) + (1 - x)flMi[4>M^{x) + </'Mi(2;)]}</'M2(y)^Ta(^^) a;, y), 

i-— Fm Fi Fa/2 / / dudxdy m%(j)M{u) 
4 Jo jo Jo 

{{u + x + y- 2)?nB0A/j(x) + (1 - x)/iA/j0^.j-^ (x) - </>|'^^ (x)]}0a,/2 (i/)/i^f,(u, y) 



(A.3) 

for non-factorizable emission contributions with the {V — A) x {V — A) effective four quark vertex, 
and 



pNMiM2 
'■LRa 



(M) 



1NM1M21 

:Ra 

ICf 



-I rl rl 



Fm Fmi Fm2 / / dudxdy 171^(1) m{u) 

4 Jo Jo Jo 

{/iMa fJ-Mi {[{u-x-y + (2;) + (n + X - y - [y) 

-[{u-x-y+ l)(t)\i^ {x) + {u + x-y- 1)(/>1/^ {x)](fii^ (y)} 
+ (u - y)mB /ijv/a [0M2 (y) - ^Ih iy)]^k'h {x)]hTa{u, X, y) 



rpNMlM2 

■^LBh 



(M) 



IC'f ^ ^ ^ 



1 /-l rl 



JO Jo 



dudxdy m\(j)M{u) 



{/iMa /^A/i { [(n - X + {x) + {u + x + y- 2)<j)lj^ {x)]<l)l,j^ (y) 

+ [{u-x + y)cPlj^{x) + (u + x + y- 2)4>lj^{x)]cl)lj^{y)} 

+ {u + y- l)mB ['/'a/2 (y) + ^M2(y)]^Mi{x)^hTbiu,x,y) 



(A.4) 



for non-factorizable emission contributions with the {V — A) x (V + A) effective four quark vertex, 
and 



rpNMlM2 I 

-^SP 

rpNM-ihh 
^SPa 



(M) = T^^^^'^^{M)+r^^l''''\M) 



NMiNh, 



(M) 



.ICf_ 



Fm Fm, Fm2 [ 
Jo 



1 fi fi 



JO 



dudxdy m\(f)M{u) 



pNMi_M2 
-SPb 



(M) 



{{u + x-y- l)mB0A/i(2;) + (1 - x)/iMi ['/'a/i (^) " (pAh(.x)]}4'M2{y)hTai^i ^^v) 
ICf ''^ ''^ 



Fm Fm, Fm2 I / I dudxdy m^cpMiu) 
4 Jo Jo Jo 



{{u + y- l)mB(/>A/i(x) + (1 - x)^a/i[</'a/i(2;) + (l)M,{x)]}4>M2{y)hTbi'^^^^y) (A.5) 



for non-factorizable emission contributions with the {S — P) x (S + P) effective four quark vertex. 
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We now present the results from annihilation diagram contributions, 
Al^f'^M) = Fm Fm, Fm, j j j dudxdy ml^u) 

{(1 - y)m%(t>M2{y)4>Ah{x) + 2fiM2 MA/i[(2 - y)4>li^{y) + y(fii2{y)WMS^)}hAa{x,y), 
AlLb^'{M) = Fm Fm, Fm^ C C C dudxdy mlMu) 

4 Jo Jo JO 

{xml(j)M2{y)(t>Mi{x) + 2^A/2 ^A/i[(l + x)(j)lj^{x) - (1 - x)4>lj^{x)]4>lj^{y)}h^i,{x,y), 
Al^^^'^M) = AI^^^'Hm)+AI^,^''\M), 
Al^M^^M) = Al^^^\M\ 

AlT'^m = Al^.^'^^iM) (A.6) 

for the factorizable annihilation contributions with the {V — A) ^ (V — A) and (V — A) ^ (V + A) 
effective four quark vertexes, and 



A^SPa^^'i^) = 4^ ^A^i I I [dudxdy mlcPMiu) 

^ ^^c Jo Jo Jo 



[(1 - y)fiM2mh(y) + 't'AhiyMMiix) + '2mi'I^M2{y)4>MS^)]hAa{x,y), 

AIp^^^{M) = i--^ Fm Fm, Fm2 I j j du dx dy m%(l)M{u){2^iM24>\i^{y)(t>hh{x) 

c i/ 

+x fXMi (t>M2 (y) [0Mi (^) - ^Mi (^)] } ^Abix, y) (A.7) 

for the factorizable annihilation contributions with the [S — P) y~ {S + P) effective four quark vertex, 
and 



A^La^^^i^) = ~hjr Fmi Fm2 I I I dudxdy m%(l)Miu)^[mb + mB{u-y)]m%4>M2{y)(t>Ahix) 

c " 

+f^MifJ'M2{[-{u -x-y + l)mB(l)lj^{x) + {-u -x + y + l)mB(l)M^{x)](t)M^{y)} 

+ [(4mf, + {u + x-y- + {u - x - y + l)mB(l)M^{x)](j)lj^{y)^hAa{u, x,y), 

ANM,M2^M-^ = Fm Fm, Fm2 f f C dudxdy m|</.MH{xm|(/>M,(y)</.Mi(x) 

+MA/iAtA/2{ - [{u + x + y- l)4f^j^{x) + {-u + x-y + l)4!l4,{x)](fii^{y) 
+ [{-u + x-y + l)<i>l,^ {x) + {u + x + y- l),^,, {x)]<Plj^ (y)} y) (A.8) 

for the non-factorizable annihilation contributions with the {V — A) x {V — A) effective four quark 
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vertex, and 



ANM.m,^^-^ = 4^ ^^'^2 /' /' /' ^^^^2;% m%ct)M{u) 

^ATMiM.^^) = -^^^FmFm, Fm, r r f dudxdy m|0M(n){x^M,K,,(x) + 0l,,(rE)]</.^,^(y) 

c ■/ 

-(1 - u - y)HM2[4^M2iy) - '&2(yM^hix)}hAbiu,x,y) (A.9) 

for the non-factorizable annihilation contributions with the {V — A) x (V — A) and (V — A) x (V + A) 
effective four quark vertexes, and 



1 Cp ( 
^SPa'^^(^) = -'i'^j^ Pm Fm^ I / / dudx dy mB<i)Miu)^[mb + {x - l)mB]rn%(j)M^{y) 

c >/ 

+f^Mim2{4>Mi{x)[{u -x-y + l)mB4>M^{x) + (-M - X + y + l)r7iB(/)X/j (y)} 

+ {[4m6 - (-U - X + y + l)mij](/)^^^ (x) - (n - x - y + l)mB(/)^^^(x)}05^2(y)|h^„('U,x,y), 

ANM,M2^^^ = Fm Fm, Fm, r r r dudxdy m%<l)M{u)[{-u - y + l)ml4>MMMx) 

+MA/i/iM2{[(^^ + X + ?/ - l)0^^^(x) - {-u + x-y + l)(t>M,{x)]4>lj^{y)} 
+ [{-u + x-y + l)(l)l^_^{x) - {u + X + y - l)(l)li^{x)](l)lj^{y)}h^i,{u,x,y) (A.IO) 

for the non-factorizable annihilation contributions with the {S — P) x (S + P) effective four quark 
vertex. 

The functions with {A = a,b) from Eqs. (A.l) to (A.9) arise from propagators of gluon and 
quark, here Y = F, N denote the factorizable and non-factorizable contributions respectively, and 
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X = T,A the emission and annihilation diagrams respectively. They have the following explicit forms: 

(— — x)m^ — fi^ + ie){xm^ — mj^ + le) 

h^biu, x) 



(— u(l — x)m'^ — /i^ + ie){—urin?Q — + ie) ' 

h^^{u,x,y) 



(— — rc)m^ — /z^ + ie)((l — x)(l — n — y)m^ — + ze) ' 



h^biu, X, y) 



(— n(l — x)'m^ — fig + ie){{l — x){y — u)m^ — + ie) ' 
(x(l - y)m| - //2 + ie)((l - y)m| - m2 + ie) ' 



(x(l — y)'in'^ — fi'^ + ie){x m?^ — + ie) 

hAaiu,x,y) 



(x(l — y)m'^^ — /i^ + ie){[y — n)(l — x)m'^^ — + ie) ' 

3;, y) = — 2 I ■ ^ 2 2~r~T- (^-^^^ 

— y)m^ — + ie)[[l — u — y)x — + le) 
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